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Abstract 

We have developed an exact method to calculate the vortex-antivortex in- 
teraction energy in the anisotropic 3D-XY model. For this calculation, dual 
transformation which is already known for the 2D XY-model was extended. 
We found an explicit form of this interaction energy as a function of the 
anisotropic ratio and the separation r between the vortex and antivortex lo- 
cated on the same layer. The form of interaction energy is lnr at the small r 
limit but is proportional to r at the opposite limit. This form of interaction 
energy is consistent with the upper bound calculation using the variational 

method by Cataudella and Minnhagen. 
PACS: 74.25.Ha,74.60.-w,75.30.Gw,75.30.Mb 
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In studying the phase fluctuation effects of high temperature layered superconductors, 
several approximations of the vortex-antivortex interaction energy have been suggested 
at the level of the highly anisotropic 3D XY-model (from now on called the layered XY- 
model) or the Lawrence-Doniach model For the bare interaction energy, Cataudella 

and Minnhagen [Ij] adopted the variational method and found the upper bound of this energy. 
According to their calculation, this interaction energy increases linearly with the separation 
between the vortex and the antivortex, which is different from the logarithmically increasing 
energy of the 2D XY-model. This is not surprising since this energy at larger separations 
should be dominated by the Josephson vortex lines connecting the vortex and antivortex 
(see Fig|). 

In passing, we must clarify what we mean by a vortex-antivortex pair on the same layer 
(Fig.^). For the layered XY-model or Lawrence-Doniach model, the vortex line cannot be 
disconnected, and should either be infinitely long or a closed loop. Thus a vortex and an- 
tivortex pair on a layer should be connected by Josephson vortex lines (Josephson strings) 
residing between the layers. Recent experimental results on high temperature superconduc- 
tors were, however, explained through the interpretation that the phase fluctuations 
are associated with the layers and that major roles in phase fluctuation effects are played 
only by the vortex and antivortex on the layer. Hence the vortex-antivortex interaction en- 
ergy is reasonably defined fUJij] as the smallest energy of the very vortex loop configuration, 
which corresponds to the shortest Josephson strings. The Josephson strings just modify the 
interaction energy of the pair at large separations. Monte Carlo simulations 0-§f performed 
to understand the phenomenon of vortex fluctuation at finite temperature also support this. 

In the present paper we develop a dual transform of the highly anisotropic 3D XY-model 
and obtain an explicit form of the vortex-antivortex interaction energy as a function of the 
anisotropic ratio and r. The dual transformation of 2D XY-model which was developed 



by Jose et al. [OJ was extended for anisotropic 3D XY-model. This result is compared 
with that of the variational calculation by Cataudella and Minnhagen ||. We also discuss 
the result in connection with the high-T c superconductors. This is under the assumption 



that the Lawrence-Doniach model is reduced to the layered XY-model. The reduction is 
possible if the amplitude fluctuation is strongly suppressed and the energy associated with 
the induced magnetic field is neglected [|]. 



I. DUAL TRANSFORMATION 



We begin with the partition function 



2- 



{d9} exp 



tjn 



{d9} = Y[d6 lthn 
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- S(i\j\n) = K\\ cos (9m,j,n - Oi,j, n ) + K \\ cos {0i,^-i, n - 9i,j,n) + K ± cos (9 itjtnH - ii3> ) (2) 

where Kn = K is the intralayer coupling constant and K± = e 2 K (e <C 1) is the interlayer 
coupling constant. In an approximate connection to the Lawrence-Doniach model, 
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where ip is the Ginzburg-Landau order parameter, £y is the Ginzburg-Landau inplane coher- 
ence length. The m (M) is the effective mass parallel (perpendicular) to the plane. After 
expanding each exponential factor in Fourier series 



JK cos < 
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where I m are modified Bessel functions, and after integrating out 9ij n , we can rewrite Z 
and S as 



Z (x^2 exp 
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- S{i\j\n) = - — \m x (i + l,i\j\n)\ 2 - — \m y {i\j + l, j\n)\ 2 - ^-^\m y (i\j\n + l, n)\ 2 , (5) 

where the primed sum denotes the constraint 



m x (i+l,i\j\n) - m x (i,i-l\j\n) + m y (i\j+l, j\n) - m y (i\j,j-l\n) 



+ m z (i\j\n+l, n) — m z (i\j\n, n — 1) = 



(6) 



for all i, j and n. 

The constraint for the summation over m can be solved by moving to the dual lattice 
of the original cubic lattice (see Fig.0) and by defining another integer field i on the dual 
lattice sites 

m x (i,i+l\j\n) = e z (i+-\j+l\ n ) - e g (i+l\j-l\ n ) - e y (i+^\j\n+^) + £ v {i+~\j\n-^) 
m y (i\j,j+l\n) = £ x {i\j +~\n+ -) - 4(«|j + ^|n-^) - + + + + 
m z {i\j\n,n+l) = e v {i + -\j\n+^) - i y {i-\\j\n + \) - lM3 + \\n+\) + £ x (i\j ~\\n+\) ■ 

In terms of the field on the dual lattice, Z is written in a matrix form 
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where the bra and ket notation is defined below 
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(For explicit form of M, see Appendix HJ.) Using the Poisson resummation rule 
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it follows that 
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where the quantity Q^iji = x,y or z) is interpreted to be the vorticity in the direction of 
/x-axes 
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Since element of the matrix M depends only on the differences between the site indices, 
it is convenient to calculate Z in the momentum space where the matrix M is diagonal. 
Then we treat only 3x3 matrix M(p x ,p y ,p z ) 

z^Y.1 W ^p f- / (dp)s(p)] , / (dp) = [ " 

{Q} J 1 J 



o (2vr) 3 



(9) 



with 

- S(p) = ~ ( £(p) | M(p) | Hp) > + 2m ( Q(p) \ £(p) ) . (10) 

Now we diagonalize the matrix M(p) (Appendix |A|) with three eigenvectors \vq), | f i ) 
and | t>2 ) and three eigenvalues u)o(= 0), cji and u 2 - We integrate over ( Vq \ £), (vi\ £) 
and ( v 2 | instead of Then, the zero eigenvalue gives a constraint to the vortex 

configuration (Q(p) \ v o ) = or 

QM3 + \\n + \) - + Qv( l + \\l\ n + \) ~ QS + \\3~ l \ n + \) 

+ Qz(i+^\j+^\n) - Q x (i+-\j+^\n-l) = 0, (11) 

i.e., the vortex line is either infinitely long or a closed loop. After Gaussian integration over 
(vi\ i) and ( f 2 I i ), we finally obtain the partition function for vortices 



{Q} 

with 



Z oc Jl'exp - J (dp)S(p) 



(12) 
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- S{p) = --(2n) 2 K(Q(p) I [ — \v 1 )(v 1 \ + — \ v, } (v 2 \ \ Q(p) ) (13) 
2 \Wi w 2 

where the primed sum reminds us of the constraint Eqn.(|TT|). 



II. RESULT 

The dual transformation discussed above provides an explicit form of interaction energy 
of the vortex-antivortex pair on a layer (Fig.|I]). Just for convenience (noting that we are 
interested in the large separation limit), we assume that the vortex and antivortex is placed 
parallel to the y-axes (Q x = 0). Let the separation between the vortex and the antivortex 
be 2r. In this case, the configuration vector (Q(p) \ is given by 

( Q(P) I = \\ e+ e Z P y ( > - e ' tPZ ^ 1 - 6 ~ lPV ) 
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If we substitute this vector into ([13]), we obtain the interaction energy in units of 4irK as 

U(2r) = tfurHe) ■ (2r) + - rfp ^f^ i^sin tp) (15) 

7r Jo sin if 

where 

F e (a) = t dt ( - 1 - v/ \ T 7? ^ . (16) 

V ; Jo \Ve 2 + l-t 2 Va 2 + 1- tWa 2 + e 2 + 1 - t 2 J V ; 

To simplify the integral form of the interaction energy of Eqn.([T5|), we used the formula 

r 

sin(2r0)/sin0 = 2 ^ cos[(2s- 1)0]. (17) 

s=l 

Note that in the case of e = 0, Eqn.(|T5|) is reduced into the interaction for the 2D XY-model, 
while even small value of e induces the term linear in r. 
The asymptotic expansion for r ^> 1 (Appendix [B]) gives 

U(2r) ~tan~ x (e) • (2r) +// c + 0(l/r), (/i c = Const.). (18) 

Cataudella and Minnhagen, in their variational calculation, took a simple approximation, 
where the phases on the layers next to the layer containing the vortex-antivortex pair are 
identically zero, and found an upper bound. Their result is expected to be larger than the 
values of our exact calculation. In our notations, the interaction energy of Cataudella and 
Minnhagen is as follows (in units of AttK) 

U var {2r) ~ • e • (2r) + Const. 

This, compared with ours U(2r) ~ [e + 0(e)] ■ (2r) + Const., is larger roughly by a factor 

7T/V2. 

For small values of r, simple asymptotic form is not available. But numerical evaluation 
shows logarithmic increase with 2r of the interaction energy (Fig.[3](b)). Plots of U(2r) for 
2r < 150 and e = 0,0.01,0.02,0.03 which include the intermediate case is shown in Fig|| 
It clearly shows crossover from logarithmic behavior for small r to linear one for large r. If 
the anisotropic ratio e increases, the region of logarithmic dependence decreases. For the 
details, we draw U-log(2r) plot at Fig.|](b). 
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The modification of the bare vortex-antivortex interaction, in consequence, leads to the 
correction to the two-dimensional Berezinskii-Kosterlitz-Thouless (BKT) transition fT7j-|20f . 
Hikami and Tsuneto |2]J studied this effect by the renormalization group (RG) analysis. In 



their study, they assumed that the two dimensional RG equation |20[ is valid, but subject 
to the cutoff 1/e because the vortex-antivortex interaction energy is logarithmic only at 
distances smaller than 1/e (in units of the lattice constant) p2| . According to their result, 
for example, the true transition temperature T c is shifted from the two dimensional transition 
point Tbkt due to the small interlayer coupling by an amount of order (7r/|lne|) 2 . This 
result has also been confirmed by the Monte Carlo simulations 



III. CONCLUSION 

We developed a dual transformation of the highly anisotropic 3D XY-model to study 
the bare vortex-antivortex interaction. We found that this dual transformation method pro- 
vides an exact form of vortex-antivortex interaction energy as a function of the anisotropic 
ratio and the separation between the vortex-antivortex pair. This form of this interac- 
tion is in good agreement with a zero temperature variational calculation by Cataudella 
and Minnhagen. The correction to the two dimensional BKT transition due to the small 
interlayer coupling was also briefly discussed. 



ACKNOWLEDGEMENT 

We wish to express appreciation for the financial support of the Korean Ministry of 
Science and Technology, and the Korean Ministry of Education. 



APPENDIX A: 

In this appendix, we summarize the informations about the matrix M in Eqn.(|3|) and 
Eqn. fllOl) . First we define some short-hand notations for convenience: 
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V(p M ) = 1 - exp(-i^), 
A(p M ) = 2(1 - cosp M ), /i = or z 
&(Px,Py) = A(p x ) + A(p y ) 
A(p x ,p y ,p z ) = A(p x ) + A(p y ) + A(p z ) 
A e (p x ,p y ,p z ) = A(p x ) + A(p„) + t 2 A{p z ). 

The Fourier transform of M^(ii'\jj'\nn') is defined as 

M^Uyflnn') = J (dp) e +^') P * e +iU-j>y e +i(n-n') Pz M ^(p x ,p y ,p z ). 

And the momentum space 3x3 matrix M(p x ,p y ,p z ) looks like 

M = e~ 2 M l + M 2 , 

' A(p„) -V(^)V*(^)0^ 

Mi = 



A(Py) -V( Px )V*(Py) 

-V*(^)V(^) A(p x ) 




Mo 



\ 



A( Pz ) -V(p x )V*(p z ) 

2A(p 2 ) -W(p y )W*(p z ) 

-V*( Px )V( Pz ) -V*(p y )V( Pz ) A(p x ,p y ) 

We need eigenvalues and eigenvectors of M(p). Below the eigensystem follows 



V 



/ 



u = 0, wi = A(p x ,p y ,p z ), u 2 = e 2 A e (p x ,p y ,p z ), 





' V(p x ) N 




^( Px )V*( Pz )) 






\vo ) oc 


V(p,) 


, | Di ) OC 


V( Py )V*(p z ) 


, | w 2 ) oc 


-V*(p x ) 




v V(p*) J 




V -A(p x ,Py) ) 




V J 



where | t>o ), | fi ) and 1 t>2 ) are to be normalized. 



APPENDIX B: 



In this appendix, we prove the asymptotic behavior in Eqn.([H|) of U{2r) for the limit 
r 3> 1. At first, we show that the function F e (a) in Eqn.(|T5D vanishes with a — * as 



F e (ot) < Aa 2 \ lna|, (A> 0). (Bl) 

And then, we show that for U(2r), the correction to the linear behavior in r is at most 
0(l/r). 

Resorting to the mean value theorem in the interval [0, a], there exists < c < a such 
that 

F e (a) = F e (0) + a ^-F e (c). 

dc 



Then the inequality follows: 



r-l 

FJoi) = ac I eft 



o 



+ 



( c 2 +1 _ t 2)3/2 ( c 2 + e 2 + X _ t 2)3/2 



rl ( 1 I 

<ac / eft - ■ + 



o " Vc 2 + l-t 2 c 2 + e 2 + l-£ 2 



- ac lvTT? ln Vi + ^-J - (0<e<o<<6) 

V2v / TTo 2 Vl + tt 2 -!/ 
~ a 2 | In a | 

Now we decompose G(a) = F e (a)/a 2 into diverging part Gca v {a) and smooth part 
G smt h{a). Then the above argument shows that the divergence of Gdi V (a) as a — > is 
at most logarithmic. In consequence, the integration 

/" 7r /2 7f 77" j 

/ dip sin 2 (2n/) llnsin^l = — In 2 • — 

Jo ' '4 8 2r 

leads to the simple asymptotic form 

mt/2 

/ efy> sin 2 (2r^)G di „(sin(^) = Const. + C(l/r), (r > 1). (B2) 
Jo 

Finally noting that 



/ dip sm 2 (2rip)G smth (smip) = Const. + 0(l/r), (r > 1) (B3) 
Jo 

completes the proof. 
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FIGURES 



FIG. 1. Vortex-anti vortex pair. In a three dimensional lattice, a vortex line is either infinitely 
long or a closed loop. For highly anisotropic 3D XY-model, however, the above configurations is 
most the important. This configuration including the Josephson string is called a vortex- antivortex 
pair. (Minnhagen and Olsson, 1991) 

FIG. 2. Dual lattice of cubic lattice. The auxiliary field I is defined at each of the face center 
in the cubic lattice. The constraint of Eqn.(6) to the field m is automatically satisfied for any value 
of the field 1. 

FIG. 3. Numerical evaluation of the interaction energy. The part (a) shows linear dependence 
on 2r for r>l. At the part (b), the crossover from logarithmic behavior for small r to linear one 
for large r is clear. For a comparison, the decoupled case (e = 0) is also plotted with filled circles. 
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